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^^ We extend and explore the general non-relativistic effective theory of dark matter (DM) direct 
i~^ detection. We describe the basic non-relativistic building blocks of operators and discuss their sym- 
O metry properties, writing down all Galilean-invariant operators up to quadratic order in momentum 
Q^ transfer arising from exchange of particles of spin 1 or less. Any DM particle theory can be trans- 
it lated into the coefficients of an effective operator and any effective operator can be simply related 
i__i to most general description of the nuclear response. We find several operators which lead to novel 
_^ nuclear responses. These responses differ significantly from the standard minimal WIMP cases in 
^ their relative coupling strengths to various elements, changing how the results from different exper- 
iments should be compared against each other. Response functions are evaluated for common DM 
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lO targets - F, Na, Ge, I, and Xe - using standard shell model techniques. We point out that each of the 



nuclear responses is familiar from past studies of semi-leptonic electroweak interactions, and thus 



^-J potentially testable in weak interaction studies. We provide tables of the full set of required matrix 

PsJ elements at finite momentum transfer for a range of common elements, making a careful and fully 

. . model-independent analysis possible. Finally, we discuss embedding non-relativistic effective theory 

. ^ operators into UV models of dark matter. 



1 Introduction and Summary of the Effective Theory 

The nature of dark matter is a fascinating mystery that continues to be unsolved. Direct 
detection experiments offer the possibility of determining the precise interactions of a dark 
matter (DM) particle with nuclei. The experiments use different targets, potentially testing 
various types of interactions between the dark matter particle and the nucleus. Previously, 
there have been attempts to exploit these differences between targets, in order to reconcile 
potential signals seen at some experiments [H E] with the absence of a signal at others. 
However, studies are typically performed in a model-driven way, with the goal of putting 
constraints on a specific particle model. Recently, [3] took a different approach, by considering 
the leading non-relativistic operators coupling DM to nuclei, and placing bounds on their 



coefficients. This was done in the context of elastic scattering. In this paper we would 
like to extend this approach by going beyond the leading operators to include the full set of 
possible operators in the non-relativistic theory, including momentum- and velocity-dependent 
operators as well. The goal of this study is to identify all possible elastic DM-nucleus response 
functions that may be exploited by experimentalists to characterize DM, and to relate these 
responses to the underlying effective theory operators that mediate the DM-ordinary matter 
interaction. 

Models with momentum-dependent interactions have received some attention lately, as 
they have helped lessen the tension between apparently conflicting direct detection experi- 
ments in El [HI [TJ E]. However, from a bottom-up point of view, their appeal is much more 
general. At the present moment, almost nothing is known about the non-gravitational inter- 
actions of dark matter with the Standard Model, and in general, assumptions about couplings 
are driven almost totally by appeals to minimalism or specific models of the electroweak scale. 
Such principles are not necessarily a good guide as to what we should expect for the nature 
of dark matter, and if the dark matter is instead as complicated as the Standard Model itself, 
then we may expect much richer possibilities for its structure and interactions. In particular, 
if the dark matter is composite, like most of the visible matter in the universe is, then one 
should expect dark matter form factors related to their compositeness scale. From this per- 
spective, momentum-dependent interactions are a compelling and well-motivated possibility, 
since they require only a small amount of structure in the dark matter sector. The usual argu- 
ment against such dependence is that it will be suppressed, since any momentum-dependent 
terms will necessarily vanish in the limit of zero momentum transfer at direct detection ex- 
periments. This however neglects the possibility that the leading, momentum-independent 
interactions can easily be suppressed or forced to vanish, leaving the momentum-dependent 
interactions as the dominant ones. 

As a simple example of this kind of model, imagine that we have a gauge boson A' that 
mixes kinetically with the photon £kin = eF^^F''^'^ [9l HOl IHl |12j. Now, take the dark matter 
X to be Majorana, so that a charge interaction with A' is forbidden, whereas the anapole 
operator xi'^I^XW'^'u is not, and can be generated when the A' gauge symmetry is broken 
J13j . In such a case, the four-fermion operator that is generated by integrating out A' is 
the interaction eg'ex'y'^'y^X^lfj.^/''^'A- In the non-relativistic limit and at zero momentum 
transfer, this is equivalent to 

which manifestly is suppressed by powers of momenta. Another, perhaps even more mundane, 
example is for the dark matter to be a Dirac fermion that is a neutral composite particle made 
up of constituents charged under the A'^ gauge force. The charged constituents will cause x to 
have a magnetic dipole moment of order their charge Q' times the radius r of the bound object, 
i.e. £dm ~ 9'Q'i^X^^''xFlu- Integrating out the A'^ again generates a four-fermion operator 
with momentum-dependence, eg'eQ'rx(y^^q^x^ltJ.^/^A- This can be seen to vanish in the 



limit of zero niomentuin, which occurs physically because at long wavelengths the interaction 
averages over the charge of the constituents and sees only a neutral object. Some of the 
earliest considerations of such scenarios are [HI [15] . 

Rather than inventing all such possible models one by one, it is more efficient to pass 
directly to an effective field theory description. Generally, such a description is the most 
natural and efficient tool to perform bottom-up, model-independent analyses. In this case, 
the appropriate effective field theory for direct detection experiments involves a set of four- 
fermion operators for the interactions of dark matter with a nucleon in the nucleus in the 
non-relativistic limit. The full set of such operators, being higher-order in the momentum, 
have also not been considered previously. Interestingly, as we will see, some of these opera- 
tors lead to novel nuclear responses, and therefore new form-factors are needed to describe 
DM interactions with the nucleus. In particular, direct detection should include not just 
spin-independent (SI) and spin- dependent (SD) interactions, but also angular-momentum de- 
pendent (LD), as well as spin and angular-momentum dependent (LSD) interactions. Under 
this new framework, the various elements used for direct detection, couple with different 
strengths, depending on their nuclear properties. It therefore becomes important to check 
whether current direct detection experiments have a "blind spot" when combined. Namely, 
whether there are any operators (or combinations of operators), which render dark matter 
less visible to the currently available targets [16]. 

Before describing the possible nuclear responses, let us first provide a quick summary of 
the non-relativistic effective theory of nucleon-DM interactions. Since we are interested in 
elastic scattering direct detection, all effective operators will be four-field operators, of the 
form 

^int = xO^^xNOnN. (2) 

The generalization to "inelastic dark matter" would involve allowing Xi^xX2i with xi ^-^d 
X2 having different masses; such a generalization should be straightforward to include, but 
we will not consider it further here. Passing to momentum space, we will take the incoming 
(outgoing) momentum of x to be p (p') and of A^ to be A; {k'). The form of possible interactions 
is constrained by several symmetries. In particular Galilean invariance imposes that the 
only combinations of momentum that may appear are those made from the two invariants 
momentum transfer q = p' — p and relative incoming velocity v = v^An — 'VN,in- Interactions 
can contain the nucleon spin S^ and, if the dark matter carries intrinsic spin (for instance, 
if it is a fermion), then S^ as well. Because the interaction must be Hermitian, it is useful to 
work with Hermitian quantities, the complete set of which is 

iq, v-^ = v + - , Sy., Sn, (3) 

where the notation v-^ is introduced because, by energy conservation, v-^ ■ q = 0. 

We will work up to second-order in the momentum exchanged between the dark matter 
particle and the nucleus. Also, we will limit ourselves to operators which arise due to exchange 



of particles of spin one or less (i.e. at most quadratic in either Sorv). In any Lorentz-invariant 
local quantum field theory, CP-violation is equivalent to T-violation, so let us first consider 
operators that respect time reversal symmetry. These operators are 

1, S^ ■ Sn, v"^, i{S^ X q) -v, iv- {Sn x q), {S^ ■ q){SN ■ q) (4) 

■y-*- ■ S^^ -y-*- ■ Smi iS^ ■ {Sn x q). 

The operators in the first line of eq. Q are parity conserving, while those of the second line 
are parity violating. In addition, there are T-violating operators: 

iSn ■ q, iS^ ■ q, (5) 

{iSN -qjiv^ ■ S^), (iS^ ■q)iv^ ■ Sn)- 

In order to determine the interaction of DM particles with the nucleus, the above oper- 
ators need to be inserted between nuclear states. Experimentally, the relevant question is 
thus what sort of nuclear responses these operators illicit when DM couples to the nucleus. 
We find that there are six basic responses corresponding to single-nucleon operators labeled 
Mj.p^n, ^'j-p,n-> '^'j;p,n^ ^J;p,n, ^'j-pni ^"j;p,n i^ ^ur discussiou of scctiou 3. Fivc of thcse re- 
sponses {Mj.p^n, ^'j;p,ni '^"j;p,ni ^J;p,ni ^",j;p,ri) ^rise in CP conserving interactions (due to the 
exchange of spin one or less), and we therefore primarily focus on this smaller set. Although a 
certain CP- violating interaction can be viable (see section 6), finding a UV-model which will 
result in the response $[/.„„ seems more challenging. In this paper we provide form factors in 
detail for some commonly used elements, however, it is useful to have a heuristic description 
for the responses. M is the standard spin-independent response. S', S" are the transverse 
and longitudinal (with respect to the momentum transfer) components of the nucleon spin 
(either p oi n). They favor elements with unpaired nucleons. A certain linear combination 
of them is the usual spin-dependent coupling. A at zero-momentum transfer measures the 
net angular-momentum of a nucleon (either p or n). This response can be an important 
contribution to the coupling of DM to elements with unpaired nucleons, occupying an orbital 
shell with non-zero angular momentum. Finally, $", at zero-momentum transfer is related to 
(L ■ S)n,p- It favors elements with large, not fully occupied, spin-partner angular-momentum 
orbitals (i.e. when orbitals j = £ ± ^ are not fully occupied). As all these responses view 
nuclei differently, a completely model independent treatment of the experiments requires data 
to be considered for each response separately (up to interference effects). 

Our paper is organized as follows. In section 2, we describe in detail the effective field 
theory, emphasizing the non-relativistic building blocks of operators and their symmetry 
properties, and demonstrate that the operators in (IfS) describe the most general low-energy 



theory given our assumptions. In section 3, we discuss the relevant nuclear physics, and in 
particular we thoroughly analyze the possible nuclear response function in a partial wave 
basis, which is the standard formalism for such physics. In section 4, we give an overview of 
the various new nuclear responses, with an emphasis on their relative strength at different 
elements. In section 5, we summarize these results in a format that can be easily read off and 



used in analyses of constraints from direct detection experiments. In particular, section 5 and 
appendix |A] contain the approximate necessary nuclear matrix elements and form factors at 
the most experimentally relevant elements. In section 6, we discuss possible models leading 
to operators which have not been considered previously. We conclude in section 7. 

2 Non-relativistic Effective Theory 

2.1 Preliminaries 

Consider the following non-standard dark matter example |1]. Let dark matter be a complex 
scalar particle that is a bound state of two fermions ipiip2-, and with compositeness scale 
A. To couple this to the Standard Model, introduce a new Uil) gauge boson A'^ of mass 
tra that kinetically mixes with the photon, C D eF' F^^ . If the fermions have equal and 
opposite charge under the A gauge field, then will be neutral. However, it will interact 
with A through the lowest- dimensional interaction that is not forbidden, which in this case 
is a charge radius coupling: g' d^ (pd" (p* F'^^ / Is? . Integrating out the A^ generates, at leading 
order in momenta, the following interaction: 

4^P[X]9'^0*0Ar7'^iV, (6) 

where in this case the nucleon A^ is the proton, and the brackets [^,y] indicate the anti- 
symmetric component. In the limit of small momentum transfer g, the v = 0,/i = 1,2,3 
terms dominate, and one is left with the interaction 



C h)irrN-N~), c="-^^^¥^. (7) 
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Here, N ,(f) are non-relativistic fields involving only creation or annihilation fields, i.e. 

This example illustrates a few points that will be useful to keep in mind when we turn to a 
systematic description of the full effective field theory. First of all, the leading interaction in 
this case is momentum-dependent, the leading standard interactions having been eliminated 
by charge assignments. Second, it is just the first term in an expansion in powers of q over 
the compositeness scale A. In order for this expansion to make sense, q must be less than A 
over the range of momentum transfers relevant at direct detection experiments. 

The minimum possible cut-off Ayy on the effective theory is dictated by the relevant 
experiments, and must be at least as large as the experimentally probed region of momentum 
transfer q between dark matter and atomic nuclei. Direct detection experiments directly 
measure the recoil energy Eji of atomic nuclei, and for a target nucleus of mass rriT this 
corresponds to a specific momentum transfer q = y/2mTEji. Consequently, the momentum 



transfer of any event is known up to any uncertainty in the recoil energy and possibly the 
identity of the atomic element. In general, then, every experiment has a maximum momentum 
transfer that it is sensitive to. This follows from the fact that larger momentum transfers 
require the dark matter to be incident with greater velocity, at least fmm = Q'/2/iT- The 
velocity distribution of the dark matter halo is expected to fall off exponentially at around 
V ~ 10~^, and essentially shut off completely at the escape velocity Wesc ~ 2 x 10~^. Since 
the reduced mass /xy is always less than the target mass my, and typical target masses are 
rriT ~ 100 GeV, momentum transfers will rarely if ever be larger than 

g^ax ~ 200MeV. (9) 

So, in order for the effective theory to be a reliable description of direct detection experiments, 
one should have Auv > few ■ gmax- Larger cut-offs are of course allowed, though they will have 
correspondingly smaller cross-sections for dark matter scattering with nuclei. 

2.2 The Effective Theory 

We will now explore the effective theory in more detail, describing the essential ingredients 
and the full set of possible non-relativistic interactions. The kinetic action is just the usual, 
non-relativistic form: 

£kin = 2m^0+(y)U^-^J0-(l/). (10) 

By momentum-conservation, the momentum transfer q is both 

q = p'—p = k — k'. (11) 

There are several important symmetries that restrict the possible form of interactions. The 
first of these is Galilean invariance, which is just a constant shift in all velocities. Thus, all 
momenta must appear through Galilean invariant combinations. Between p, k and q, there 
are therefore only two independent momenta that can arise in any interaction. It is easy to 
see that the momentum transfer g*is Galilean invariant, as is the relative incoming velocity 

V = V^,in - VN,m, (12) 

which is just the velocity of the incoming dark matter particle in the nucleon rest frame. The 
final kinematic constraint is energy conservation. This is easiest to impose by passing to the 
center-of-mass system, where the total kinetic energy is E = ^f^NV^ei^ with /i^v = m'^+m ^^^ 
dark-matter- nucleon reduced mass and v^c\ = v^ — vn- For the initial state energy, this is 
just -Ein = ^fiNv"^, whereas for the final state energy it is Eout = ^fJ'N^v + -2-)^. Imposing 
i?in = -Eout is therefore equivalent to 

v-q = --^. (13) 

2hn 



The next major constraint is Hermiticity of the interaction. This is essentially equivalent to 
crossing symmetry, because Hermitian conjugation exchanges incoming for outgoing particles, 
i.e. (0~)^ = (f)'^ ■ Consequently, the momentum transfer g'is effectively anti- Hermitian, and it 
will be more convenient to work with the Hermitian operator iq. Under exchange of incoming 

t, r. ^ ^ , q 



and outgoing particles, v does not have definite parity: v — )■ ■y'x.out — "WAf.out 
we can easily construct a similar quantity that is Hermitian: 



v^ 



ii-N 



V^ 



q 



2/i 



N 



However, 



(14) 



The reason for this notation is that, by the energy-conservation condition above, v^ ■ q = 0. 
Finally, we must include the particle spins. In the relativistic limit, this is just the familiar fact 
that four-fermion operators can contain 7 matrices. In the non-relativistic limit, we can write 
down the dark matter and nuclear spins S^ and Sn as operators directly. Different possibilities 
for the spin of the dark matter are thus treated in a unified way. If dark matter is a spin-1/2 
particle, then these spins operators are simply ^a, where a* are Pauli sigma matrices, acting 
on the X aiid N spinors; for vector dark matter, they are spin-1 representations of the angular 
momentum generators J* acting on the x vector; and for scalars, they simply do not appear. 
These are invariant under Hermitian conjugation, so we have for our complete set of Galilean, 
Hermitian invariants the following: 



iq, 



s. 



XI 



Sn- 



(15) 



In addition to the above symmetries, there are strong constraints on violations of CP sym- 
metry. Since ultimately our non-relativistic theory must be embedded in a Lorentz invariant 
quantum field theory, this is equivalent to T symmetry. Spins behave like angular momen- 
tum, and thus change sign under T. Also, velocities all change direction under T, so v-^ and q 
change sign as well. Finally, although we will not impose P as a symmetry, it will be helpful 
to classify all operators according to whether they are even or odd under P. In this case, spins 
do not change sign, whereas v-^ and q do. Thus, the complete set of Galilean invariants has 
the following transformation table: 





t 


T 


P 


s 


+1 


-1 


+ 1 


iq 


+1 


+ 1 


-1 


v^ 


+1 


-1 


-1 



Since we are interested in elastic scattering direct detection, all effective operators will be 
four-field operators, of the form 



jC^ 



int 



X+O^X-N+OnN~ = O x^X-N+N' 



(16) 



Furthermore, the momentum-transfer-squared g^ is a completely invariant scalar quantity 
that depends only on dark matter kinematic quantities, and thus if O is an operator allowed 



by all symmetries of the theory, then q^^O is as well. It is therefore natural to classify all 
such operators as a single one with a g^-dependent coefficient, or form factor: 

coO + c^g^O + C4g^O + . . . = FoU^O. (17) 

Massless mediators can be incorporated by including a Fq ~ q~^ term, though strictly speak- 
ing this is not an effective operator. A related point is that at the upper range of momentum 
at experiments, the pion should be included in the effective theory and x-x-tt couplings al- 
lowed. For instance, if the underlying DM model contains couplings such as XltJil^xJ^ ^^ 
DM to the axial current J^ = ig'^^'j^r^q, then the effective theory will couple x's to pions 
due to the overlap of J^^ with vr. Such interactions would contribute to dark matter-nucleon 
scattering through t-channel pion exchange at tree-level, effectively producing Fq oc 21^2 
form factors in x-^-A^-A^ interactions. 

So far, we have mainly discussed momentum scales. In addition, there is an energy scale 
associated with the scattering process, of size Uq ~ q^ /2mT < 200 keV. This is usually negli- 
gible, as the binding energy uj of nucleons is about 10 MeV per nucleon for most elements, and 
inelastic transitions are kinematically suppressed. However, for nuclei with small splittings 
~ ujq between the ground state and an excited state, it could affect direct detection rates. 

We are now ready to present the possible non-relativistic interactions. The general La- 
grangian is 

^'-t = E E^S^'^^.X+X^iV+iV-, (18) 

N=n,p i 

with the following set of operators. Of the T-even operators, we have 

1. P-even, S'^^-independent 

Oi = l, 02 = {v^f, 03=iSN-{qxv^), (19) 

2. P-even, S*^- dependent 

0, = S^-Sr,, 0, = tS^-iqxv^), 0, = {S^-q){SM-q), (20) 

3. P-odd, S*^- independent 

Oj = Sm-v^, (21) 

4. P-odd, S';^-dependent 

Os = S^-v^, 0, = iS^-{SNxq) (22) 



In addition, we also have T-violating operators: 
5. P-odd, S"^- independent: 



6. P-odd, S'y-dependent 



0,0 = ^SN ■ q, (23) 



On = iS^ ■ q. (24) 



It is convenient to separate the operators as we have done above because each of these six 
groups of operators will not interfere with each other. In addition, there are four operators 
that are products of the ones above: 

O10O5, C?io08, 01103, and On07. (25) 

With these, the above operators provide the most general effective theory at the dark matter- 
nucleon level that can arise from exchange of a spin-0 or spin-1. In the completely general 
effective theory for elastic scattering, one would relax this condition and include arbitrary 
powers of v and 5*^, which would allow products of the operators we have written here and 
one additional operator O12 = S^ ■ {Sjy x v-^). For instance, OyOg is a local operator that we 
have not written down above. However, quadratic powers of Sn and beyond (and S*^ as well, if 
X is spin- 1/2) can always be reduced to at most linear powers by using the multiplication table 
of sigma matrices. In appendix [C| we present the non-relativistic reduction of all relativistic 
operators arising from a spin-0 or spin-1 exchange (or more precisely, with at most a single- 
index field exchange at tree-level) in terms of the local interactions above. The product 



operators in eq. (25) are seen to arise from a spin-1 particle coupling to fermion bilinears of 
the form Nd^'y^N, which, for model-building concerns to be discussed in section [6| we will 
not focus on further. Nevertheless, it should be kept in mind that the general effective theory 
possible without any such restrictions contains these operators. 

In order to obtain the size of scattering cross-sections relevant to experiments, we need to 
evaluate matrix elements of the nucleon-level operators from the effective theory inside of a 
target nucleus. From the point of view of the effective field theory we have constructed, an 
atomic nucleus is a heavy, many-body bound state of nucleons. For the purpose of computing 
nucleon matrix elements inside such a nucleus, it is important to separate out v-^ into a term 
vj^ that acts on the coherent center-of-mass velocity of the atomic nucleus as a whole, and a 
term vj^ that acts only on the relative distances of the nucleons within the nucleus. We can 
write 

I 

{Vx,m + ^^x.out - VN;m " ^Ar,out) = ^^T + ^7V, (26) 



where 



2 



7,{Vx,m + ^^x,out - VT,m - VT,out) = Vt + (27) 



acts only on the center of mass motion of the nucleus (here, Vt = v^,m — VT,m is the incoming 
dark matter velocity in the lab frame). Also, vj^ is just 

vif = --{vNAn + VN,ont), (28) 

but where the "uVs act only on the separation distance between the nucleons. The reason for 
this separation is that vj^ and vj^ behave qualitatively differently, and have parametrically 
different sizes. The former is determined by the kinematics of the dark-matter-nucleus scat- 
tering process, and does not require any detailed knowledge of the internal structure of the 
nucleus. Its approximate magnitude is given by 

rriTVT ~ q- (29) 

Indeed, for elastic scattering, by kinematics, vt must be strictly greater than ^^, and the 
event rate tends to be suppressed by the halo distribution if vt is significantly greater. As 
stated above, the typical size for vt is ~ 10"'^. 

On the other hand, v^ depends on the internal distribution of nucleons in the nucleons 
and thus is determined by 

rriNvjj ~ q. (30) 

This will lead to a relative kinematic enhancement of mT/mN = A for vj;^ compared to vj=^. 
In many cases this is cancelled by the fact that vj^ tends to sum coherently over nucleons 
whereas vj^ often does not, making the two terms comparable, but we will see some important 
exceptions. 

At low momentum-transfer, the internal structure of atomic nuclei can be summarized in 
just a finite number of macroscopic quantities. In the case of the standard spin-independent 
interaction Oi or spin-dependent interaction 04, these are the atomic number A and charge 
Z or nucleon spin expectation values (S'„), (Sp), respectively. However, there are many more 
possible macroscopic quantities that appear associated with our full table of interactions than 
just these usual ones. Furthermore, at finite momentum-transfer, there are multiple possible 
form factors associated with the nuclear responses that are required for calculating event 
rates. Still, there are fewer independent nuclear responses than the full set of operators in the 
effective theory, so that a small number of plots can roughly capture the full range of possible 
models. 

In order to obtain these nuclear responses, one needs detailed input from nuclear physics 
on the wavefunctions of nucleons inside the nucleus. In the next section, we will therefore turn 
to the standard framework for the results of these computations, and a thorough discussion 
of possible nuclear responses. We will provide a mapping of the operators in the effective 
theory onto the nuclear response functions. 

In section |4| we will discuss how these nuclear responses favor different elements. The 
reader whose immediate goal is to use the resulting form factors in order to compute specific 
experimental event rates will find the relevant results summarized in section [5] and appendix 

10 



A. 3 A series of nuclear physics calculations of moderate complexity have been carried out in 
order to illustrate the kinds of variations among nuclear target responses one should expect, 
given the unknown nature of DM-nuclear interactions and the range of effective theory possi- 
bilities. While we would characterize our nuclear structure calculations as reasonable - based 
on the shell model, using realistic interactions that have been "vetted" in related electroweak 
studies, and employing bases of reasonable size (ranging up to ~ 0.7M Slater determinants, 
after applications of symmetries) - we also hope our results will motivate others to train even 
more sophisticated nuclear structure technology on this problem. For example, the relevant 
isotopes of Ge span a region in neutron number where a sharp spherical-to-deformed transi- 
tion occurs, accompanied by fascinating quantum-mechanical level-crossing phenomena and 
associated sharp changes in proton and neutron spectroscopic factors. This paper will pro- 
vide those specialists with tools necessary to tackle such problems an important additional 
motivation for undertaking new work, its relevance to ongoing experiments that address one 
of the most important open questions in particle astrophysics, the nature of DM. 

3 The Nuclear Responses 

The non-relativistic effective theory treatment of dark matter responses connects naturally to 
the standard language of multipole expansions for nuclear electroweak responses that we sum- 
marize here. Such expansions allow one to exploit nuclear selection rules based on rotational 
invariance, parity and time reversal. We specialize here to the case of elastic dark matter 
interactions, as the energy transfers in dark-matter scattering generally preclude inelastic 
excitations. (However, the extension of the multipole formalism to inelastic dark matter in- 
teractions is straightforward and will be presented elsewhere.) The good approximate parity 
and CP of the nuclear ground state then impose important new selection rules on the possible 
elastic operators, restricting the multi-polarities that contribute as well as their interference. 

3.1 Nuclear Charges and Currents 

The usual construction of coordinate-space nuclear charge and current densities in electroweak 
interactions begins with a covariant interaction that is reduced to produce the needed non- 
relativistic operators. The effective theory approach significantly simplifies this analysis, and 
also provides important guidance to those who might want to follow a model-dependent 
analysis based on some specified covariant interaction. Our starting point is the interaction 

C-EFT = ctil + 02^"*" ■ v'^ + usSn " (q X v'^) + a^S^ ■ Sn + io-bS^ ■ {q x v'^) + a&S^ ■ qSN ■ q 

-l-ayS'iv ■ y-^ + agS'^ ■ ■w"'" + iagSy^ ■ (Sn x ^ + io-ioSN ■ q + io-uSx ■ q (31) 

As we have discussed previously, the Hermitian velocity f -*- can be divided into a target center- 
of-mass piece vj^ and components vj^ associated with the relative velocities of target nucleons 
(and thus with the A-1 Jacobi momenta). In combination with nuclear spins, these two 
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velocities generate interactions that are separately invariant under Galilean transformations. 



Before any models of nuclear charges and currents are introduced, Eq. (31) shows that from 
the available nuclear degrees of freedom - nuclear spins and relative momenta vj^ - one can 
construct the nuclear "charges" 1, vj^ -vj^, and Sn ■ vjj- that transform under parity and time- 
reversal as even-even, even-even, and odd-even, and nuclear "currents" vj^, Sn, and Sn x vj^ 
that transform as odd-odd, even-odd, and odd-even. Given our assumption of a nuclear 
ground state with good parity and CP, this leads us to conclude that there must be six 
independent nuclear response functions corresponding to the even multipoles of a vector-like 
charge operator Mjm, the odd multipoles of axial longitudinal L^jm, axial transverse electric 
Tf^, and vector magnetic T^^ operators, and the even multipoles of vector-like longitudinal 
LjM and transverse electric Tj\.^ operators. To go further - to provide explicit forms for these 
multipole operators - we must make nuclear model assumptions. Our construction defines 
the spins and momenta Sn and vj^ as the local operators associated with nucleons. This 
one-body definition is the most common starting point for nuclear physics calculations. 

We do an explicit example - the axial-charge operator O-j = Sn • v-*- - to demonstrate 
the procedure for separating f-*- into its vj^ and vj^ components and for constructing the 
associated nuclear operators. First, in the elementary-particle point limit for a nucleus - 
where the nucleus is characterized only by its macroscopic quantum numbers of charge, spin, 
and isospin - we have 

1 ^ 1 

I ~* I ~* I -L X — ^ I -L 

O7 = v -Sn > ^j.-Sn = Vj.--y a{j) where Vj. = - (^7^,^ + i^x.out - ^r,in - ^>r,out) • (32) 

point / ' ' Z 

(A factor of one-half is introduced in relating the nuclear spin Sn to nucleon Pauli spin 
operators.) The center-of-mass nuclear velocities VT,in and ■ur.out can be obtained by averaging 
over the velocities of the nucleons, 



-^ A 1 "^ 

VT,in = -T^VN,m{j) and VT,ont = -T^VN,ont{j) ■ (33) 

3=1 i=i 

The nuclear model-building assumption is that the underlying nuclear charge and current 
operators are one-body and local, in this case the sum over the individual (symmetrized and 
thus Hermitian) axial charge operators. This operator can be explicitly separated into its 
center-of-mass (which contributes to ■u^) and intrinsic components (the vj^ contribution) 



j=l j=l j>k=l 

A [A ^. .. 



2^T 



j=l 






(34) 

intrinsic 



The first (target recoil) term is already identified in Eq. (32) while the second provides an 



explicit definition for the contribution associated with the A-1 relative Jacobi three-momenta 
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and thus with vj^. While in principle the first, explicitly Galilean invariant form for the 
intrinsic axial-charge operator - a two-body operator in relative coordinates - could be used 
in calculations, in many cases the simpler one-body form can be employed provided the center- 
of-mass motion of the nucleus is properly treated. Thus the interpretation of the intrinsic 
subscript on the axial-charge operator above is an instruction that such steps should be 
taken, if this form of the operator is employed. This is typically done by working in a 
translationally-separable (e.g., full shell) harmonic oscillator Slater determinant basis, then 
numerically forcing the center-of-mass to reside in the Is state: the limitations of this approach 
are discussed in the appendix. 

An advantage of the effective theory treatment is that it immediately identifies the trans- 
lationally invariant recoil axial charge contribution to dark matter scattering proportional to 
vj^. For elastic scattering, this is the only contribution of the axial-charge: matrix elements 
of the intrinsic operator vanish for even multipoles by parity and for odd multipoles by time 
reversal. This result would be more difficult to obtain in conventional treatments that begin 
with a covariant interaction. Such calculations would need to extract the recoil term from 
the axial-charge operator (in contrast to having it manifestly in the point-nucleus limit of the 
effective operator), a task often requiring the combining of charge and current contributions. 
For example, consider the example of a V-A four-fermion contact operator between dark mat- 
ter and a nucleus, ip^'j^ipx 4'N'j'^l^i'N- Defining l^ = ipy^'j^ip^, the charge and three-current 
contributions to scattering are 






Y.-ij) 



PfU)+PiU) , VT,in + VT.^ 



2m7v 



-> 



E 



a(i) + 






PfU)+PiU) 



2ni 



N 



intrinsic 



tV 



-I ■ i^Niibi^N -^ 



^xM + ^■ 



X,out 



E^(^)' 



(35) 



where we have inserted the spin operator as the non-relativistic limit of the axial three-current 
operator. Indeed we get the right answer: summing the two terms yields a contribution pro- 
portional to vj^ as well as the intrinsic operator. But in contrast to the effective theory 
treatment -where the target contribution is immediate from the point-nucleus limit and ap- 
pears as one term - a certain degree of care is needed to locate and regroup terms into Galilean 
invariants. 



We can now handle the general case of Eq. (31), first arranging the various terms as 
follows 



Let = lol + lo [-24 • Sn] + h ■ [25^] + hi ■ [-v 

'Pi + Pf 



ki + i-^""^"' 



2m 



N 



a + k ■ a + I 



M 



2m 



N 



Ni+h- [2i 4 ^ Sn 
■ Pi+Pf 



+ 1 



E ■ 



2m 



X a 



N 



(36) 



where the coefficients of the charge (/q), axial charge (Iq), axial vector (Z5), vector magnetic 



(/jv/), and vector electric [Ie) densities, determined from Eq. (31), will be given below. This 



expression follows Eq. (31) exactly apart from one simplification, the elimination of the 
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term proportional to v^ ■ vj^ within 02- This interaction transforms as a parity- and time- 
reversal-even charge of o{v/c)'^ ~ 1%, and thus will be overwhelmed by the o{v/c)'^ coherent 
spin-independent response, if the latter is present. Furthermore, from a model-building point 
of view, it is difficult to see how one could generate the former while avoiding the latter, 
without significant fine-tuning. 

The charge and current operators can be transformed to coordinate space via the substi- 
tution 



Pi+Pf _^ 
2m]sf 2m 



— ( — VS(x-Xi) + S(x-Xi)-v] 
UN \ I % J 



Thus we determine the Hamiltonian density 

A A 

Uet{x) 



^/o(»i) 5(f -fi) + ^^c 



° ^'^ 2M 

A A 

+ y^ h{i) ■ S{i)S{x - Xj) + y^ ImJJ) ■ 

i=l i=l 

1 



-i^< 



(7{i)5{x — Xi) + 5{x — Xi)(7{i) 



(37) 



i^. 



2M 



-^r5{ 



X — Xi) + 5{x — X 



M 






E[^ 



2M I 



V j X (j{i)5{x — Xi) + 5{x — Xi)a{i) x Vj 



(38) 



where the dark-matter amplitudes /o(i) and l{i) appear within the sum over nucleons because 



we will allow the various couplings in Eq. (31) to have a nontrivial isospin dependence, e.g. 



ai — )■ (aj + a\T^{i)) (so that a\ = a\ = a\j2 will correspond to a coupling only to protons of 
strength a\ while a^ = — cti = ai/2 will correspond a similar coupling only to neutrons). 



The Hamiltonian for Eq. (38) has the familiar form 



dx e"*"-^' 



UJM^\P{X)\J^M,) - I ■ {J,MS{x)\J,M,) 



(39) 



where q is the three-momentum transferred from the nucleus to the scattered DM particle. 



One substitutes Eq. ( 38 ) into Eq. ( 39 ) and uses the spherical harmonic and vector spherical 
harmonic identities 



^iq-Xi 



5^v^[J]z^jj(gx,)r;o(fi.J 
I 
oo ^ 

5^V^ [J] t'-^^jj{qXi)Yj^{n^:, 



j=o 

/ oo 



exe 



iq-Xi 



J=0 






;J-2 



Ajj(gx,)rJVi(fi.,) + y X jj{qx,)Y}j,{Sl^,) 



X 


= 




(40) 


, A = 


±1 



to project out charge multipoles and longitudinal, transverse magnetic, and transverse electric 
current multipoles, respectively. This defines the operators that generate the nuclear form 
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factors describing the scattering of dark matter from nuclei. Here [J] = -\/2J + 1 and ex, X = 
1, 0, —1, are spherical unit vectors defined with respect to a quantization z-axis along q = q/q. 
As detailed in the appendix, the multipole operators transform simply under parity and 
time reversal, allowing one to exploit selection rules to simplify the diagonal nuclear matrix 
elements of interest, assuming CP- and parity- violating components in nuclear ground-state 
wave functions are negligible. Consequently, averaging over initial nuclear spins and summing 
over final, one finds the general form of the dark-matter elastic scattering probability: 



2Ji + 



^ J2 UMf\H\JM)\' 



An 



Mi,Mf 



2J, + 1 



Y, UW k ■ q ^'M U)\' 



.J=l,3„... 



oo ^ 

+ V |(J,||/oMj(g)||J,)p + |(J,||4-g^<|."(g)||J,)p 



J=0,2,... 

+ 2Re 

,2 °^ 



J,\\ lE-q—^"{q) mm k Mj{q) ||J,)=* 



+ ^ E ((^di4<i'',(g)iiJ.)-(j.ii4<i'S(g)ii^.r-i(^.ii4-g$j(?)ii^.)r) 

A'" 7_o A 



J=2,4,... 



+ Y. {^{^J^W^M ^j{q) ||J.)-(J.||rM Aj(g) ||J,)*-|(J,|| lM-q^j{q) m? 



J=l,3,... 



'Af 



+ 2 (<^^ll ^^ ^'-^(^^ ll"^^^ • ^"^'11 ^^ ^'^(^) ll"^^)* - 1^"^^" ^' • ^ ^'^(^) ll"^^)!') 



+ 2Re 



iq-{-m hi —Aj{q) \\J,) X {JiW k ^'Aq) WJiY 
rriN 



(41) 



All nuclear matrix elements are intrinsic: contributions proportional to vj^ reside entirely in 
the dark-matter amplitudes /q, h, Ie, and Im, by virtue of the Galilean invariant effective 



theory. In Eq. (41) || denotes a nuclear matrix element reduced in angular momentum. The 



expression is somewhat schematic in that 



[J,\\ I Oj{q) U) = (J.ll 5^/(^) Oj{qx,) U), 



(42) 



i=l 



The notation is a reminder that the dark matter amplitude in general cannot be moved out- 
side the nuclear matrix element because that amplitude may contain several effective theory 
couplings with different isospin dependences. (There are many cases where this expression 
does factor, however, and we give the simpler form appropriate for those cases below.) 



Equation (41) shows that there are six distinct nuclear response functions governing dark- 



matter responses, corresponding to six single-particle operators. Each of these operators is 
familiar from standard treatments of weak interactions [171 HH] or, in the case of $' and $", 
from extensions [19] that have been made to account for currents of order 1/m^. They are 
constructed from the Bessel spherical harmonics and vector spherical harmonics, MjM{qx) = 
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Jj{qx)YjM{n,) and Mji = jL{qx)Y^ 



JLM, 



MjM{qx) 

q 



AjMiqx) ^ Mfj{qx)--V 



^Miiqx) = |^VMjM(gx)|-a=[J]-i{vOTTMJ5+i(gx) + v^MJS_i(gf)}-<T 
$';^(gf) = « (-VMjMiqx)] ■ (a x -V^ (43) 



q / \ q 

The multipole operators have been defined to have a simple behavior under time reversal, 
transforming with a ±1, as discussed in the appendix. Time reversal and parity impose im- 
portant constraints on allowed responses: for reasons noted previously, there is no elastic 
dark matter coupling to the intrinsic axial charge density (though the axial charge contri- 
bution due to vj^ remains and contributes through the spin density). Those responses that 
do appear involve sums over either even or odd multipoles, again because of the parity /time 
reversal constraints. The long-wavelength limits of these operators, showing explicitly the 
character of the nuclear response (charge or current; transverse magnetic, transverse electric, 
or longitudinal; vector-like or axial- vector like) are given in Table 11} 

Table [l] includes the standard spin-independent response governed by even multipoles of 
the (generalized in isospin) charge operator Mjm, and two spin-dependent responses that, 
though proportional in the long- wavelength limit, are characterized by different nuclear form 
factors and couple to dark matter in distinct ways, with Sjj^,^ being longitudinal and T,'jj^ 
transverse. The associated form factors involve sums over all allowed odd multipoles. There 
is a third J = 1 response, a transverse magnetic response governed by odd multipoles of 
AjM- This and two other new responses are explicitly associated with nuclear substructure. 
AjM is generated by the convection current (the nucleon velocity term), which in the long- 
wavelength limit produces a coupling to the nuclear orbital angular momentum operator 
£{i). This provides a third interaction - an interaction like the two spin-dependent ones - 
that will transform under rotations as (JjI Ja/I Jj). The two other responses arising from the 
constituent nature of the nucleus transform as longitudinal and electric projections of the 
density ~ 6{x — Xi)a{i) x V. The first of these is quite interesting as its long-wavelength 
limit produces a scalar proportional to the spin-orbit interaction (t{i) ■ i{i) as well as a tensor 
contribution. The full form factor involves a sum over all even multipoles of $". Because of 
the leading a(i) ■ i{i) contribution, this response is present for all nuclei (that is, regardless of 
ground-state spin, like the usual spin-independent charge coupling), but because of its form 
factor (leading-order behavior is proportional to g/rrijv) and spin-orbit nature, its properties 
are quite different from those of the usual spin- independent scalar operator Mqq. $gQ can be 
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Response x 



4-!T 

2Ji+l 
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Long-wavelength 
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Type 
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oo 
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oo 
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rriN 
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J=2,4,. 
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^AiMiqXi 
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■iriN 
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mjv 



"^00 (?^^; 
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2V37r 
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(^iM{i) 



0-lM[i 



2mAr\'6n 



'imUJ 



3mivv47r 



o- U 



m,jvv SOtt 



mjvv 207r 



[a;^ (a(i) x fV)i]2M 

[Xi ® ((T(i) X }V)i]2M 



MjM ■ Charge 



^JM • Axial 
Longitudinal 

-'■JM '■ Axial 
Transverse Electric 

ry-imag 
-'■JM '■ 

Transverse Magnetic 

LjM '■ 
Longitudinal 



-'-JM ■ 

Transverse Electric 



Table 1: The response dark-matter nuclear response functions, their leading order behavior, 
and the response type. The notation ® denotes a spherical tensor product, while x is the 
conventional cross product. 
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important in heavy nuclei because a{i) ■ £{i) produces a coherent isoscalar contribution over 
closed spin-orbit partner shells, e.g., the closed I/7/2 shell for Ge isotopes and the closed lgg/2 
shell for Xe or I. The operator $ " transforms as the longitudinal projection of a vector current 
(denoted Ljm in Table [I|. Just as in the case of the nuclear spin density, where two distinct 
spin- dependent form factors are generated, corresponding to the axial-like longitudinal and 
transverse electric nuclear responses (denoted Lj and Tf^ in Tabled]), the density responsible 
for $ J also generates a transverse-electric response (denoted TJ^^v/) that consequently transform 
as a J = 2 operator in the long-wavelength limit. From the specific form, [xi^{a{i) x ^V)i]2m, 
one can see that this operator is closely related to $qo — )■ <j{i) ■ i{i) and is in fact proportional 
in the long- wavelength limit to its tensor partner ^2m- Thus the relationship of $j to $j - 
Lj and Tf operators - is analogous to that of S'j to Sj - Lj and Tf^ operators - except that 
the transverse nature of ^'jm excludes the possibility of a J=0 multipole. We will see below 
that $jAf is an exotic response, arising only for dark matter with unusual couplings. 

The three composite operators arise from the fact that there are interesting current den- 
sities in the nucleus that can mediate dark matter interactions, but fail to have the proper 
parity and time-reversal properties to contribute in the point-nucleus limit. As momentum 
transfers in dark matter interactions are not small compared to the inverse nuclear size, these 
new responses can be numerically quite important. More important, we will see below that 
in many cases, these new responses can provide the dominant coupling of dark matter to 
nuclei, depending of the effective theory operator. The long-wavelength limits of these three 
new responses are determined by operators that transform properly under parity and time 
reversal because a factor of q- Xi has been convolved with the underlying bare nuclear densi- 
ties. Consequently these operators have a leading-order form-factor behavior proportional to 
q/rriN and an explicit dependence on Xj, and thus on the nuclear size. 

3.2 The Effective Theory Content 



Thus we can proceed to the dark-matter physics, which is encoded in the amplitudes /q, I 



Im and Ie that are determined by our effective theory through Eq. (31 ). We find 



^0 = («? + 0^x3(2)) -i(gx S;^) -{T^ (a^ + 4r3(i)) + 4c- ^^T («8 + 4^3(^)) 
+ iq- S^ {a\^ + a\^T-i{i)) 

h = ^ iqxvr («3 + 4^3(0) + S^ (0-4 + alr-iii)) + 5^^ ■ g g (0° + 4^3(0)" 

+ v^ (a° + a\T3{{)) + iqx S^ (a" + 09X3(2)) + iq {a^Q + al^r^ii)) 

Im = iqy< S^ (ttg + 05X3(2)) - S^ (ag + al)T3{i) 

Ie = ^g(a° + 4r3(z)) (44) 

We observe that there is no coupling to the Tf nuclear density associated with the tensor 
operator $j=2- ^3 generates a nonzero Ie, but it is longitudinal. None of our eleven effective 
theory operators generates a transverse component to Ie- This point relates to our decision 
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to limit the effective tlieory to interactions tliat could arise from spin-0 or spin-1 exchanges. 
Indeed, we pointed out the existence of an additional operator linear in dark-matter and 
nuclear spins 

012 = Sx ■ (Sn X v^) (45) 

that would need to be included if we were to relax conditions on the nature of the exchange. 
With the inclusion of this term, 



Ir -^ 



q {al + 4x3(2)) + iS^x(«i2 + «12^3(0) 



(46) 



leading to a contribution to the sixth potential nuclear response, the tensor one, governed by 
$'. O12 would also produce a recoil contribution to the spin-dependent response functions. 
Another example of a term that produces a contribution to $' is the four-fermion tensor/axial 
tensor interaction 



-"tensor 



-ipx^fiuiar - ^«at75)V'x ^Ncr^^iaT - iO'ATl5)'^N 



.N 



(47) 



Inspecting Eq. (44) one sees that there are several common situations in which our general 



result for the scattering probability, Eq. (41), can be simplified by factoring the dark- matter 



amplitudes from the nuclear matrix elements. If we are interested in any one interaction Oi, 
then clearly its associated isospin dependence can be written aj(l + ajr3). The overall strength 



Cj could be associated with the dark-matter amplitudes in Eq. (44), while the isospin factor 



could be included in the definition of the single-particle operators of Eq. (43). Alternatively, 



several couplings might be nonzero, but all might share a common behavior in isospin, e.g., all 
interactions coupling to protons. The overall couplings could again be incorporated into Eq. 



(44), with the isospin dependence (1 + Ts{i))/2 absorbed into the single particle operators. In 



such cases the scattering probability simplifies, taking the form 



+ 
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+ 2Re 



tq-ifM X 5) {J^\\ -^Aj \\J^){J^\\ T.'j{q) ||J,)* | 



(48) 



As the inultipole expansion is conventionally done in a coordinate system aligned along 
g*, another useful result is the expression for the Hamiltonian in a rotationally invariant form, 
in the long wavelength limit. One finds 



i=l 

-i^ilM X q) ■ I{i) + ^(4 ® q)2 ■ \xi ® [a{i) x \V{i)]^ ] (49) 



J 2 



The third and fifth terms, proportional to /s and {lE®q)2-, can each be divided into longitudinal 
and transverse electric pieces, which would be associated with distinct nuclear form factors 
once one goes beyond the long wavelength limit. 

3.3 Response Function Evaluation 

The various response functions described above were evaluated in the shell model for several 
of the key isotopes now used in dark-matter detectors. Calculations were performed for ^^F, 
23Na, 70,72,73,74,76Qg^ i27i^ ^^^ I28,i29,i30,i3i,i32,i34,i36xe. Response functions were evaluated by 

summing over the contributing isotopes, weighted according to their natural abundances. 
(For F, Na, and I, there is a single stable isotope.) Consequently, while all isotopes take part 
in scalar or spin-independent responses, only those with ground-state spins >l/2 (^^F(l/2'^), 
23Na(3/2+), ^3Ge(9/2+), ^^Tj (5/2+), ^^^Xe{l/2+), and i3iXe(3/2+)) contribute to J = 1 
(or spin-dependent) responses, and only those with spins > 1 (^^Na, ^^Ge, ^^^I, and ^^^Xe) 
contribute to the J = 2 tensor response. By defining the nuclear responses per target atom 
for a detector made up of unenriched isotopes, we take into account the reduced efficiency of 
detectors for J = 1 and J = 2 responses due to noncontributing isotopes. The Ge calculations 
include the five stable isotopes, while the Xe calculations summed over the seven principal 
isotopes, ignoring the trace (< 0.1%) contributions from ^^^'^^®Xe. 

As our focus is a broad survey - to understand the degree to which targets can vary in the 
relative sensitivity to dark matter, given the broad range of response functions that may gov- 
ern that sensitivity - the structure calculations we undertook were limited to relatively small 
bases, and thus should be considered exploratory. They were performed in m-scheme bases 
on which we placed a limit of no more than 0.65 million Slater determinants (after application 
of symmetries like time reversal to reduce basis dimensions). The sd-shell calculations for ^^F 
and ^^Na are then unrestricted. The interaction used was that of Brown and Wildenthal [20] . 
The Ge isotopes were treated in the standard 1/5/22^1/22^3/2 1(79/2 model space above a ^^Ni 
core. The basis truncation was based on limiting occupation of the lgg/2 shell to no more 
than two nucleons above the minimum occupation for all isotopes. An interaction developed 
by the Madrid/Strasbourg group was used [211 ES] . These Ge model spaces take into account 
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some of the polarization effects that accompany the rather sharp spherical-to-deformed tran- 
sition that occurs near neutron numbers 40-42 (^^Ge, ^^Ge), though we recognize the need for 
further expansion of the basis in follow-up calculations. Some of the motivating physics for 
more ambitious calculations - particularly the rather complex interactions among spherical 
and deformed 0+ bands in the even isotopes - is discussed in |23]. 

The ^^^I and Xe isotopes were treated in the 851/22^3/22^5/2 1(77/2 l^ii/2 model space above 
a ^°°Sn core. The interaction used was one developed some time ago by Baldridge and Vary 
[21] and employed in double beta decay studies in this mass region. The interaction is based 
on a G-matrix from the Reid soft-core potential augmented by phenomenological pairing and 
multipole forces. While the ^^^Xe and ^^^Xe calculations were unrestricted, significant trun- 
cations became necessary for lighter Xe isotopes where the neutron occupation of the l/iii/2 
shell drops. The bases for 128,130,132-^g ^^^ i27j ^gj^g limited by fixing the l/iii/2 occupation 
to the minimum allowed nucleon number. Basis for the odd-neutron isotopes ^^^Xe and ^^^Xe 
were further restricted by limiting valence protons to the energetically favored 2^5/2 and 15^7/2 
shells, and by requiring neutrons to fully occupying these same shells (a choice that preserves 
good isospin). Less restrictive calculations can and should be done, but are beyond the scope 
of the current survey. 

From these wave functions the ground-state to ground-state one-body density matrices 
can be generated 



1 



[J][T] 



'^t^M,m = n^i^{J^^Tf\\ a„|®ci|/,| _:;j.;T,) (50) 



''I I 



J J:T 



where a\i3\^mj,mt = (~1)'''^ mj+1/2 '"ta|/3| _mj,-mt, \oi\ denotes are nonmagnetic quantum num- 
bers, ® denotes a spherical tensor product, and :: indicates reduction in both angular mo- 
mentum and isospin. As we include all contributing multipoles, density matrices are needed 
for < J < 2 Jj and for T = 0, 1. They provide the single-particle amplitudes needed for 
evaluating many-body matrix elements of any one-body operator, 

A 

i=l \a\,m 

= V2j2^'r'(}»im{\a\\\d, 

A 

{jf,T, ■■:■. ^6j(^)r(z) ;; j,-t,) = Yl ^ir='(l«U/3|) (bl ^^Ojt\\ 

(51) 

where {\a\ \\0j\\ \(3\) is a single-particle space/spin matrix elements reduced in angular mo- 
mentum. By adopting a harmonic oscillator single-particle basis, these matrix elements can be 
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evaluated analytically, yielding forms ~ e~yp{y) where p{y) is a polynomial in y = (g6/2)^, h is 
the oscillator parameter and q the magnitude of the three-momentum transfer. Consequently 
analytic expressions for nuclear form factors can be provided, so that response functions can 
be easily evaluated for changing experimental conditions, such as different WIMP masses, or 
for other choices of h. Our numerical results were generated with the choices 6=1.833, 1.835, 
2.108, 2.282, and 2.292 fm for ^^F, ^SNa, the Ge isotopes, ^"^^ and the Xe isotopes, respec- 
tively. Harmonic oscillator matrix elements for four of the operators arising in dark-matter 
elastic scattering, MjM{gx), AjM^qx), Sj^j(gif), and T,jj^{qx) can be evaluated from a pub- 
licly available Mathematica script [25] (as these operators also arise in standard treatments of 
weak interactions). A generalization of this script that includes the two additional operators 
^'jj^{qx) and $j^(gx) is available from the authors. 

4 Comparing the novel responses for different elements 
of interest 

The most important lesson that the general EFT of dark matter-nucleon elastic scattering has 
to teach us is exactly what are all the phenomenological properties that distinguish nuclear 
recoil rates at different experiments. Specific models or effective operators for dark matter 
interactions will be proportional to some particular combination of these properties, which 
can enhance or diminish the relative sensitivity of different experiments. Coupling through 
atomic number A or charge Z in the standard spin-dependent case or through the proton 
or neutron spin in the standard spin-dependent case are by far the best-known examples. 
However, as we have seen in the previous section there are other possible nuclear responses. 
In this section we would like to explore the less familiar responses, A and $". 

4.1 Ap and A„ 

The A responses in the zero momentum transfer limit simply measure the nucleon angular mo- 
mentum content of the nucleus. Therefore, elements which have an unpaired nucleon (either 
n or p), in a non s-shell orbital are favored. For the proton response, this includes ^^Na and 
i27j i9p^ whose proton is approximately in the 2si/2 orbital, is disfavored, however. Among 
odd-neutron nuclei of interest, ^^Ge and ^'^^Xe exhibit the strongest response, with other iso- 
topes less favored. In Fig. [l] at finite momentum transfer, g, we provide a comparison of the 
strength of the response for various elements (integrated over a representative range of g, and 
weighted by the natural abundances of isotopes). As explained in section 2, the A response 
receives a kinematic enhancement of A that is competitive with the coherent enhancement 
factor associated with the center of mass of motion of the nucleus. Thus, the A response 
contribution to operators such O5 and Og can be important, and can become dominant for 
elements with unpaired nucleons in large angular-momentum orbitals. To illustrate this point 
we have included in Fig. [2] a comparison of the A response to the standard SI response, M, 
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as they occur for operator Os- Finally, we would also like to point out, in case spectral data 
becomes available in future experiments, that as a function of recoil energy (or momentum 
transfer), the A responses has different behavior from either the more standard SI (M) or the 
SD (S', S") responses. This is shown in Fig. [s] 

4.2 $;' and $'^ 

The $" responses at zero momentum transfer are sensitive to a product of the nucleon spin 
and its angular momentum. The dominantly coherent part of this product, is the scalar inner 
product {L-S)n. For completely filled angular momentum orbitals, this dot product vanishes. 
Namely, when all 2(£ + 1) states of the spin-aligned {j = i + ^) subshell and all 2i states of 
the spin-anti-aligned {]=(- — \) subshell are occupied, this dot product vanishes. In general, 
however, as the C^\ orbitals have different energies and so the highest occupied orbital for 
a given element will not be filled. Let n±{t) be the approximate occupation numbers of the 
£ ± ^ orbitals. In terms of these, the dot product is proportional to (£ + 1) n+(£) — £ n_(£). 
Usually, the least energetic orbital will be filled first, and so one expects a mismatch between 
n±{€) of order I (the strong spin-orbit nuclear force in nuclei moves the spin-aligned orbit 
lower in energy). Consequently, ((L- S)) ~ ^highest ^^^ most elements. The $" responses tend 
to favor heavier elements, as these have larger i orbitals not fully occupied. Much as in the 
case of the delta responses, $" receives a kinematic enhancement of A, and can be important. 
It is the dominant response for the operator O^. As we will see in section |6j models which 
contain O3 can also typically contain :^Oi. Due to the ilighest enhancement, for heavier 
elements, the $" response can easily be of order the M response in a large portion of the 
parameter space of such models. In Fig. [T]we show the $" responses for the various elements. 
This response is particularly interesting in the context of light DM, as the sodium coupling 
strength can be more than ten times bigger than that of fluorine (see Figjll). 

5 Presentation of Results 

Ultimately, we are interested in the prediction for the differential scattering rate -^^ (per 
unit time per unit recoil energy) with respect to nuclear recoil energy. This is related to the 
differential cross-section through 

dR _ / p^rriT da ^ 



dEji \ jsj^m^v d cos 9 

where p^ is the dark matter density, and (...) indicates average over the halo velocity distri- 
bution!^ The differential cross-section depends on the matrix-elements-squared in the usual 

^ This formula follows straightforwardly from the recoil energy in terms of the velocity and scattering angle 

2 

in the center-of-mass frame, Eji = ^^v'^{l — cos 6'), and the fact that the rate per unit time is i? = (n^av). 
A canonical review of dark matter direct detection is 26 . 
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Figure 1: Size of integrated form factors J^ ° \(ldqF{q^) for different nuclear responses, 
weighted by the natural abundances of isotopes. 
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way, 



da 



dcos9 



1 1^1 



\M\' 



2]^ + 1 2j + 1 ^-^ 327r (m^ + 771^)2' 



(53) 



spins 



where we have averaged over 2j^ + 1 and 2j + 1 initial dark matter and nuclear spins, and 
summed over the final spins. The matrix-elements-squared in general contain interference 
terms between the different operators, and this leads to a large number of possible different 
form factors. A general Lagrangian of the form 
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Y,ct^Ot^ + cf^O\ 



(54) 



i=l 



will therefore lead to a matrix-elements-squared that can be written 






iN)(N')rj,{N,N'),2 2 



iv'.q'). 



(55) 



spins 



i,j=l N,N'=p,n 



where the form factors F}^ ' (g^) are defined as the coefficients of the q's in this relation. 



u 



and are defined to be symmetric in {i,N) ^ {j,N'). We give approximations for them at 
the most relevant nuclei in appendix |X} We have factored out the generic kinematic term 
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to the standard spin-independent (Fm) ones, for protons in iodine and fluorine. Ratios have 
been normahzed to 1 at y = 0. 
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(integrated against the velocity distribution) for different nuclear responses. 
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-^ which arises due to the conventional relativistic normahzation of states. Because the 
operators fall into sectors that do not interfere with each other due to symmetry and WIMP 
spin, only a few of the off-diagonal (z 7^ j) form factors -Fj/s are non-zero. In summary, the 
master formula for the detector event rate -^-^ (per unit time per unit detector mass per 

unit recoil energy) in terms of the form factors F^- ' and operators coefficients c\ in the 
effective theory is 

dRp _ T^r Px^T I 1 V^ v;^ IN)(N')j^(,N,N') 



Nr 



\X E 4"^"'^ir^"'\-^a), (56) 



where Nt is the number of target nuclei per detector mass. 

Finally, let us give the explicit connection between these general form factors and the 
convention for the form factors in the standard spin-dependent case, which is our 04. In 
the standard spin- dependent interaction, the usual convention is to write the coefficients of 
the operators Of and O^^ in terms of isospin-respecting and isospin-violating parameters 
oo = o„ -|- ap,ai = ttp — an respectively, which are related to the coefficients cf ,cf in eq. 



(I54|)by 

C4 = (32\/2)mAfm^Gi?aAr, {N = n,p) (57) 

where Gp is Fermi's constant and DM spin j^ = ^ is assumed. It is also conventional to 
define form factors Sqq, Su, Sqi. 

spins spins spins 

(5J 
Consequently, they are related to our F44 ' s with DM spin j^ = ^ according to 

-^44' = ^7^— — jt('S'oo + 5*11 + 5*01), F44' = ————(S00 + Su- Soi), (59) 

fS'^^ = ^i4"^ = ^(^^(5oo-5n). (60) 

6 Models 

One of the more interesting operators we have found is O^, as it leads to non-trivial dependence 
on the nucleon angular momentum. Since this operator is somewhat unusual, let us provide 
a sketch of a model where such an operator might arise. O3 can be obtained from the non- 
relativistic limit of (see appendix) 

xYx^if^txuq^'N -)■ (2m^)g^ + Smj^m-^iv ■ {q x Sm) = 2m^q'^Oi - SniNm^Os. (61) 

Our goal will then be to provide a model which generates this type of operator and describe 
its parameter space. An example of a model of this type is one which contains a new Dirac 
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pair of colored fermions, U, If and D,D, which also carry charges iQn and iQd under a new 
gauge boson, A'^. We can imagine that there is some UV sector which couples the left and 
right handed quarks of the SM {q and u,d) to our new fermions. If this UV sector respects a 
new-particle parity symmetry, then upon integrating the UV sector out we get a Lagrangian 
containing the terms: 

'^ = j^[yu{Uu){Uu)^ + y,{Dd){Ddy] (62) 



+ 



l^Kiq^Huplf'U^ + c.c) + y'.iq^H^dpb'^D^ + c.c)] + ■■■ 



A3 



Here a,/3 are two-component, left-handed, spin indices, and we assumed some structure in 
the UV theory which couples t/'s to the up sector and D's to the down sector. |j 

Upon integrating out the heavy fermions, and including the Higgs vev, w, the induced 
coupling of quarks to the new gauge boson take the form 

C = fiij^^ry-2Fl,^d''[yuQuUR'j''uR + ydQddR7''dR] (63) 

+ K2j^^F'^^[MuQuUa^''u + MdQdda''''d], 

where Ki are numbers of order one, and the Mi are of order the heavy fermion masses. 
We assume that the dark matter is also charged under the new gauge boson, and so after 
integrating the gauge boson out, there is an induced local interaction between dark matter 
and nucleons of the form 



^ = (4/)2A^M2^ ^7.X (ci, d'ijn^p) + ci„ d'in^'^n) + ■ ■ ■ ) (64) 

+ ^^]^X1,X d. (C2p pa^'p + C2. na^-'n) + ■■■, 

where the various c's depend on the parameters of the model and can be adjusted separately, 
for example by changing the fermion masses and charges. Adjusting the masses and charges 
(both of which are radiatively stable) one can reduce the part of the interaction which couples 
to q'^Oi, and vary the coupling to protons and neutrons through O3 independently. Note, 
that for the heavier elements the $" response can be of the same order as the M response 
(see Fig. Isl). Therefore, in a significant portion of parameter space of this model, the new 
LS-response can dominate the more standard SI response for some elements. 

Another interesting operator is the CP violating spin operator Oio = iq ■ Sn as it couples 
to only a portion of the nucleon spin (resulting in the response E"). A model generating such 



^For example, the four-fcrmion terms above can arise from the exchange of scalar fields charged under 
hyper- charge and A' . 



29 



1.2x10" 
1.0x10" 
8.0x10"' 
6.0x10"' 
4.0x10"' 
2.0x10"' 



q^M„ 



,17 1.3x10" 



lqX"l 



5.0x10" 1.7 xlO'-* 



Na 



?'<&;' 





6.4x10" 




6x10"* 








5x10"* 








4x10" 


|4qV'l 




3.8x10" 




3x10" 






2x10" 






i 




1x10" 


9.6x10"^ 




^ 




9.9x10'" 8.9x10'- 












F Na 


Ge 


I 


Xe 





cfM„ 



2.0x10"- 
1.5x10"- 
1.0x10"- 
5.0x10"- 



2.0x10" 



2.1x10" 



4x10"' 
3x10"' 
2x10"* 
1x10"* 



lq'M„"| 



5.1x10" 



6.1x10'-' 2.0x10'* 



Na Ge 

|4(?Vl 



9.4x10'-'' 



6.4x10" 1.5x10'^ 



2.7x10'-'' 



Figure 5: Comparison of the $" form factor to the standard spin-independent {Fm) form 
factor, with relative coefficients as given by the moment operator of the model in section [6j 



an operator is as follows. Consider a scalar whose couplings violate CP, but preserve CP 
in any given sector. After EW breaking, its couplings take the form (schematically): 



Upon integrating cj) out, we get the following Lagrangian 



£ 



VxVq- ■- 
mi 



{A-Kfml 



xxm + /« 



(47r)4m2 



qq iq-f^q + 



(65) 



(66) 



VxVq 



The first term will lead in the non-relativistic limit to Oiq and can dominate the interaction 
between x and quarks for small Yukawa couplings, while the last term produces CP violating 
effects in the SM. In particular, the last term can mediate direct CP violating decays in 
the kaon system. Experimental constraints place a bound on the mass and couplings of 
> (300 GeV)^. This still allows for significant rates at direct detection experiments. 
Finally, let us discuss some of the restrictions we made on the effective theory. Specifically, 
we have neglected operators that require contracting two indices from the dark matter part 
of the operator with the nucleon part. In order for such interactions to be generated by tree- 
level exchange of a mediator, the exchanged field would have to be spin-2 or higher (or, at 
least a non-standard anti-symmetric B^'^ spin-1 field). Higher-spin fields are subject to strong 
theoretical constraints and typically arise as composites fields with size comparable to their 
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mass; this is the case for instance with spin-2 resonances of QCD. There is nothing wrong with 
such particles, however they will usually be accompanied by spin-0 or spin-1 resonances as 
well, and it is unlikely that the higher-spin exchanges will dominate the interactions. We have 



made one additional truncation which is to neglect the product operators in eq. (25). This is 
motivated by the fact that any spin-1 exchange leading to such operators would necessarily 
couple to the dark matter sector as a CP-even field and to the nucleons as a CP-odd field, 
or vice versa. In a theory where CP is broken, this is not necessarily forbidden. However, 
it requires more work to see if such a framework can be UV-completed, in particular in a 
manner where the spin-1 mediator couples dominantly to a current in the UV in CP- violating 
way. 

7 Discussion and Future Directions 

The exact nature of dark matter remains as yet unknown, and little can be said for certain 
about its interactions with the Standard Model. Especially when comparing the results from 
different experiments, it is therefore prudent to keep an open mind about what form such 
interactions can take, and to avoid prejudices about underlying models when possible. The 
most efficient tool for separating out assumptions about UV physics and parametrizing the 
relevant low-energy possibilities is effective field theory, which we have applied here to the 
direct detection of dark matter. Effective field theory has been considered in the context, 
the most thorough analysis probably being [3]; however, the full set of possible interactions, 
including all derivative couplings and momentum-suppressed interactions, as well as the pos- 
sibility of interference between different operators, has not previously been explored. Since 
direct detection experiments are sensitive to interactions with finite momentum transfer, it 
is entirely possible and well-motivated for momentum-independent operators to be absent, 
and for momentum-dependent interactions to be responsible for the leading source of direct 
detection scattering. It turns out that the full set of possibilities is much richer than the 
standard cases, and can favor atomic nuclei in a qualitatively different way. 

We have systematically constructed the low-energy, non-relativistic effective field theory 
describing direct detection scattering. We have paid special attention to the basic non- 
relativistic building blocks, that connect directly to relevant experimental observables, and 
to their symmetry properties. This simplifies the identification of all possible non-relativistic 
operators and makes their interpretation more physical, as well as explaining patterns in the 
form of non-relativistic operators that arise when taking the non-relativistic limit of standard 
four-fermion operators. It furthermore eliminates the need for embedding the field operators 
in full representations of the Lorentz group, and thus gives a unified description of all possible 
spins for the dark matter particle. 

This effective field theory is necessarily for interactions between dark matter and nucle- 
ons. To make contact with experiment, one requires the matrix elements of these operators 
between atomic nuclei. This marriage of an effective theory treatment of DM interactions 
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with a treatment of the nuclear response shows that there are six independent nuclear re- 
sponse functions characterizing DM elastic interactions with nuclei, and that these response 
functions are associated with six single-particle operators having the requisite transformation 
properties under parity and time reversal. The new response functions are associated with the 
nuclear convection current and related spin-velocity currents that depend explicitly on nu- 
clear compositeness. In addition, two interference terms arise, adding additional complexity. 
Such complexity is helpful, providing more diagnostic handles for experimentalist to exploit, 
as they seek to determine the nature of DM. Our effective theory of DM interactions that 
include CP preserving exchanges of spin-1 or less utilizes five of the six possible responses. 

We have implemented this formalism by completing shell-model calculations of moderate 
complexity for several of the critical targets- ^^F, 23Na, 70,72,73,74,76Ge, ^^n^ I28,i29,i30,i3i,i32,i34,i36xe 
- using realistic effective interactions. Operator matrix elements were evaluated in a harmonic 
oscillator basis, a choice that allows one to express the needed form factors as polynomials 
in the square of the three-momentum transfer. The results show a wide range of sensitivities 
to underlying effective theory interactions, highly dependent on the choice of nuclear target. 
In a follow-up paper, we will analyze experimental constraints on the full effective theory 
using these form factors, though as more accurate form factors for heavier elements become 
available in the literature, such analyses should become increasingly reliable. 

Aside from the improvement of nuclear form factors, the results here can be continued 
in several directions. Models with inelastically-scattering dark matter can be well-motivated 
and lead to qualitatively interesting predictions, and it would be useful to extend the effective 
theory analysis to include such operators and any additional possible nuclear responses. Also, 
while some combinations of the operators here arise easily from UV models, for others it is 
less obvious whether or not tuning is required. It would be interesting to understand better 
if natural models for the full effective theory can be constructed or not. 
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Appendix 

A Nuclear Matrix Elements and Form Factors 

A.l Partial Wave Decomposition of Operators 

To make connection with results in nuclear physics, in section |3] we decomposed effective the- 
ory interactions into partial waves. Such a decomposition allows one to package the matrix 
elements-squared into reduced matrix elements through the Wigner-Eckart theorem. In this 
appendix, we will review some of the details of these computations. We will then in section 
A.2|give the form factors for the general effective theory in terms of a smaller set of indepen- 



dent form factors. Finally, in section A. 3 we will give numerical approximations for those 
independent form factors. 

We will begin with a few examples, starting with Oi = 1. Expanding e*"^'^ in partial waves, 
the corresponding DM-nucleus scattering matrix element takes the form 

oo „ 

Mrr',ss' = — 5ss' V V4^(-z)-'[J] (r | / d^xjj{qx)Yjo{q ■ x)pp{x)\r'), (67) 

where [J] = \/2J + 1. The expansion of e^'^'^ in this way is what leads to the introduction of 
the operator MjuigX^) = ij{qXi)YjM{q ■ Xi): 

{r\MjM{qXi)\r') = ld^Xi^*{xi)^i{xi)jj{qXi)YjM{q-Xi), (68) 



where ipi{xi) is just the harmonic oscillator wavefunction for the nucleon corresponding to i. 
Rotational invariance has been used in order to pick a fiducial direction for q, so that Mjm 
depends only on qx. 

Consider for the next simple example the T-violating operator Oiq = iS^ ■ q. The decom- 
position into partial waves takes the form 



oo 

e^'-^Oio = "2^ ■ V.e^"-" " "3 E ^^i^YW " V,MjM{qx), (69) 

j=o 

where we have integrated by parts. More generally, in the partial wave analysis, one treats 
vectors like S^ by decomposing A = X]a=o ±1 ^^^-^^ ^^^ using the identities in eq. (40). 



Then, an arbitrary operator of the form i ■ S]\j with £ a constant vector can be decomposed as 

f ■ S^e^^- ^IYI (v^(0'[^]) U^^JoiQS) - E ^ (ASjA(gf) + ^S',,(gf)) j , (70) 

where T,jM{qx) = MjjM{qx)-a. T^jm has the wrong parity to contribute to elastic scattering, 
so we may discard it in the following. Consequently, C4 decomposes as 

e^'^-O, - \S^-Y1 {:^^yy]) U^'jom - E ^(-^^^A + ^^'^a) 1 • (71) 
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Such manipulations can be performed for all the operators in the effective theory, for which 
one obtains the following matching: 



Oi = 


1 


EjCjMjo 


o^ = 


iSN ■ {qx v) 


T.JCJ (i$'io - (^^-x v^) ■ (Ea=±i 2-^e;(zS',,))) 


o,= 


Sx ■ Sn 


'A ■ Ejcj (eo^S% - EA=±i2-^e;(^S',,)) 


o,= 


iS^ ■ {qx v) 


(4 X ^9") ■ Ej cj (v^Mj,{qx) - ^ Ea=±i 2-ie;(AA,,)) 


o,= 


iSx-q)iSN-q) 


(?"'■ 4)1 ■ Ej cj (eo^S'Jo - Ex=±i 2-^e;(^ES,)) 


o,= 


Sn-v^ 


Ejcj {-Iv^ ■Ex=±i^-''en^'jx)) 


Os = 


Sx-^^ 


Sx ■ Ejcj [vj^Mjo - ^ EA=±i2-^e;(AA,,) 


o,= 


iSx ■ {Sn X q) 


-Wq^ Sx)-EjCj (EA=±i2-^e;(zSS,)) 


Oio-- 


= iSN-q 


-lEjCjq^'jo 


On-- 


= iSx-q 


i^Sx ■ q) Y.J cjMjo 



Here, cj = \/47i{i)^[J], and in all cases, the scattering amplitude Airr',ss' is given by ^ 
times the matrix element (s, r| . . . |s', r') of the nuclear response function in the table above. 
Our convention for the reduced matrix elements (j||Tj||j') is 



{j'm'\TjM\jm) 



V ^"^ 



f J J 

-m' M m 



{j'WTjWj). 



(72) 



For elastic scattering, Sja, CIjq and A'j^ do not contribute due to their parity, and all form 
factors can be written in terms of a small set of independent form factors that depend only 
on the nuclear responses: 



X 



{e) - 



2j 



4vr 'P^ 



1^ 
j=o 
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for X = M, S', S", A, $", are required for the diagonal matrix elements. Additionally, 
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for (X, Y) = (M, $") and (S', A) appear when there is interference between different re- 
sponses. We emphasize that these independent nuclear form factors depend only on the 
nuclear physics and are not all special cases of the form factors Fij for the effective theory 
operator coefficients q. In particular, the -Fj/s can depend on the dark matter spin. 

If one prefers a basis of isoscalar c*^°^ = c*^""-* + c'-^-' and isovector c^^^ = c^^-* — c*-"^ couplings, 
rather than the basis of neutron {N = n) and proton [N = p) couplings we have chosen here. 



then one can use an isoscalar-isovector form of the general event rate formula eq. (56), 
dRn _ P-^rrir I 1 



dEn 



N' 



327rm|m^ 



ij a,b=0,l 



(75) 



34 



with 






^;? 



4 v^J 






:i{n,n) 



+ F, 



ip,p) 



■^{p,n) 



■^{n,p) 



4 y' «j ' «j «j «j 

_ /'__p("'") _|_ _p(P'P) _|_ _p(P'") _ _p(".P) 



«J 



«J 



«J 



*J 



(76) 



A. 2 Form Factors for the General Effective Theory 
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The full set of form factors F^- ' (defined by eq. (55)) necessary for a general model can be 



written in terms of the basic independent nuclear form factors (defined by eq. (73) and (74)) 
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where C{j^) = (4j^(j^ + l)/3) is a prefactor that depends on the DM spin j^ and has been 
normahzed to C(|) = 1. All interference terms not given explicitly above can be seen to 
vanish. We have not included O2 here since it does not appear at leading order from any 
relativistic interaction without cancellations. 



A. 3 Approximate Form Factors 



Here, we provide form factors for the basic responses defined in eq. (73) and (74). They de- 
pend on the momentum transfer through the dimensionless variable y 
harmonic oscillator parameter, b 



(g6/2)^, where b is the 
v/41.467/(45A-i/3 _ 25A-2/3) fm. For the interference- 



type responses, Fj^y 7^ F^y \ so these are given separately, whereas for the non-interference- 
type ones, we have Fj 
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B Bound State Scattering Amplitudes: Momentum Space 

The calculation of scattering matrix elements for dark matter off of atomic nuclei starting with 
the nucleon-level Lagrangian is in general a complicated and difficult calculation, involving 
contributions from many different spin configurations of a multi-body system. However, in 
practice, at the end of the day most amplitudes depend on the complicated nucleus state 
through only a handful of quantities, such as the total atomic number, charge distribution, 
average nucleon spin, etc. The calculation of matrix elements for scattering of atomic nuclei 
in terms of these macroscopic quantities is a fairly straightforward application of quantum 
mechanics of bound states, which we now discuss. In the body of the paper, we have worked 
mostly with a partial wave basis of operators and a nuclear shell mode basis of states, which 
has several advantages. However, there are some advantages of a momentum-space basis for 
the nuclei as bound states of nucleons, not least of which is the fact that such a basis is more 
familiar in particle physics. Thus, it is somewhat simpler to make connection with standard 
formula for field operators and matrix elements. Furthermore, it is typically easier to separate 
out the center-of-mass motion of the nucleus in such a basis, which can make couplings to 
the coherent motion of nucleus easier to identify. 

Our starting point is the expression for the target non-relativistic bound state at rest: 

/d^k 1 

s,s- 

For simplicity, we are considering the bound state to be made of two particles, of mass irti 
and m2, and the total target mass uit is approximately mi + m2, up to a relatively small 
binding energy. The generalization to more particles is straightforward and in fact not very 
different from focusing on one constituent at a time and grouping the remaining constituents 
into a second much bigger constituent. The dependence on and sum over spins s,s' of the 
constituents will be left implicit from now on. The expression for the bound state \T{kout)) 
with non- vanishing momentum fcout is obtained simply by shifting the velocity by f = kont/fnT'- 

/dP'kd^r 1 

-— — e^^Xr) , ^ \k + rmv,-k + m2v), (88) 

where we are using the same symbol if) for the wavefunction and its Fourier transform. 

Scattering matrix elements can be computed by evaluating matrix elements of interactions 
on the states \T{k')) constructed above. Our primary concern will be to understand contri- 
butions to V from the internal motions of the nucleus, so we will consider the example where 
we are calculating matrix elements of the operator v-^ N~ N^ (y) . n Since we are dealing with 
non-relativistic physics, we have separated out the field operator N{y) into its positive- and 



■^ In order for the interaction to satisfy Galilean invariance, it is important that w^ in the Lagrangian is an 
operator that acts non-trivially on the fields N^ , so that it gets contributions not only from the nucleus center- 
of-mass velocities Uin, Wout, but also from the internal momenta fc, k' of the nucleon states \k + miv), \ — k+m2v) 
within the nucleus. Explicitly, v-^N^N+ = ^j^i^^^^^ ^ dN-N+). 
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negative-frequency parts. They satisfy [N~{y),a^j^{k)] = 0,[N~{y),a]\f{k)] = — JmT ^^^^'^ ■> 
and similarly for A^"*" = A^"''". Let us take particle 2 to be N , and particle 1 to be some other 
particle. Then, 

2mi2m]\f 

k' + k 

x( (t;i„ + v^^^))2mi5^^\{k + rriiv^^t) - {k' + miv-,^)). (89) 

niN 

The (5-function inside the integrand imposes the relation k — k' = rriiivi^ — fout), so 

k-k' + mj^{vin - t^out) = (mi + mN){v;,, - v^ut) ~ rnrivn, - v^^t) ~ hn - Knt = q, (90) 

where the approximation mi + m^r ~ ttlt is valid as long as the binding energy of the bound 
state is small. Then, making the following change of variables, 

k = Uk+^^)^ k'=kCk-^q), r-r' = x, ^=X, (91) 

the scattering matrix element reduces to 

{T{h^)\v^N-N+iy)\T{k,^,)) = e^"'^ / -d'kd'xd^Xe'''-^^'''^-'' 

rriN J 8 

X (rriTiv,^ + v,^,)r{X - ^)^(X + ^) - ^^^*(X - ^)2^x^(X + ^) 



(92) 



By integration by parts, we have traded a k for a derivative with respect to x, which is then an 
asymmetric derivative dx = \{dx — dx) with respect to X acting on the wavefunctions. Doing 
the d^k and d^x integration is now trivial. We can further simplify the expression by making 
the approximation mjv ^ f^i ~ f^Ti which is accurate to about 0{m]\f/mT) ~ 0{1/Z) for 
our physical cases of interest. Defining the density p{X) = iIj*{X)iIj{X), we have 

{T{k,,)\v^N-N+{y)\T{k^^,)) = (93) 



D^q-y 



rrirp 



mT(t;i„ + t;out) / rf3xe*'?-^p(X) -2^ / d^Xe'''-''r{X)2dx^l^{X) , 



The first term is recognized as just the bound state momentum times a form factor 



mT{v,^ + v,^,)jd'Xe'^-''p{X) = KrFiq' 



(94) 



which reproduces the well-known fact that the form factor in the spin-independent case is 
the Fourier transform of the density. The second term is more complicated, and can be 
evaluated in an expansion in small momentum transfer q. Let us define the probability current 
J = ip* {X)2dxip{X) . This is conserved diJi = 0, and as a consequence J* = dj{X'^J^) is a 
total derivative and therefore its spatial average vanishes. So, at g = 0, we get no contribution 
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from J*. The leading non-vanishing term in the small q expansion is ~ J d^Xiq^X'^JK First, 
note that X'^V-'-' = ^dk{X^X^J'') is a total derivative, which must have vanishing spatial 
average, so J d^XX'^J^ is anti-symmetric in its indices. Consequently, 

"^^ q' f d^XX'J^ = t^^q'e'^'' [ d^X^*{X){Xxpf^{X)=t^^q'L%e''K (95) 



niN J rriN J rriN 

We thus have for the full matrix element that 

{T{hr,)\v^^N-N^y)\T{K^,)) = ^^ fm^« + vl,jF{q') + ^!!^e'^''LU'S{q')] , (96) 

rriN \ rriN J 

where S{q^) is some form factor encoding the remaining subleading q^ dependence. The 
coefficient ^^^ of the final term is the effective A enhancement factor discussed in the body 
of the paper. In the relevant case of atomic nuclei, where many A^ particles are present in 
the bound state, this amplitude should be summed over all such particles. This results in an 
additional coherence factor for the first (spin-independent) term, and sends Ljv — > (-^Af) for 
the second (spin-dependent) term. 

C Non-relativistic Matching 

Let us reduce the full set of relativistic operators into non-relativistic form. For this purpose, 
we will label the incoming(outgoing) momentum of the dark matter x to be p^{p'^), and the 
incoming (outgoing) momentum of the nucleon A^ to be k^{k'^). The momentum transfer q is 
defined to be q^ = p'^ — p^ = A;'^ — k^. We will also define pt^ = pf^ -[- p't^ and Rf^ = k^ + k'^. 

The velocity v = v^in — VNin = — — = —^ — I- ^ ;r^ is defined so that it is a 

kinematic invariant {fiN is the reduced dark-matter-nucleon mass). Since there are 2 scalar 
operators {xX ^^^ Xl^x) ^^^ 4 vector operators {P^XXi P^Xl^Xi X^^^^^lvX^ ^"^^ Xl'^l^x)^ 
there are 2^ + 4^ = 20 possible combinations of operators: 
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16m%m^08 


xY-i^xNia^o.q^'N 




Wm^mNiS^ ■ {q x Sn) 


Wm^mNOg 


xYl^xNYl^N 




— IdniNm^S^ ■ Sn 


— IGmNm^O^ 


iXl^l^xK^'N-i^N 




— IGrriy^mNV-^ ■ S^iq ■ Sn 


— IGrriy^mNOioOs 


iP'^Xl^xK^.NN 




l&m%m^iq ■ S^^ 


lQm%my^On 


iP^Xl^xNia^^aq^N 


(4m^)(i 


q ■ S^){-q^ + ArriNiv^ ■ (g x ^^v)) 


Am^Oii{-q^ - AttinO^) 


tP^^Xl'xNi^^l'N 




-lQm^mN{iq ■ Sy^)v-^ ■ Sn 


—l^m-^niNOiiO-j 


Pf'Xl^xK^.N-f^N 




— IQiTLNm-^q ■ S^q ■ Sn 


— IQmNm^OQ 



Many combinations have not been included since they are equivalent by the equations of 
motion, and the above terms tend to give simpler non-relativistic pieces. The most commonly 
used such combination is the vector interaction j^, which can be written in terms of the above 
by using the Gordon identity: 

1 



N^N 



2m 



N 



-N (K^ + taf'^q^) N. 



(97) 



Note that every non-relativistic operator occurs in the above table, except for O2, which 
appear if there are cancellations in the leading pieces, for instance through the linear combi- 
nation {AmNm^xx^N ~ P'^XXK^NN). 

D Single-particle Operators 

Four of the six operators introduced in section 3 are familiar from standard treatments of 
semi-leptonic electroweak interactions [TTtlT^. MiMJqx), AjMiqx), T,'jj^{qx), and Ejj^^(gx). 
The matrix elements of these operators between single-particle harmonic oscillator states, the 
most common basis for nuclear physics calculations, can be evaluated analytically, yielding 
explicit forms for the nuclear form factors governing DM scattering. A Mathematica script 
|25j and tables [18] are available. The remaining two operators are a symmetrized form of ^'jm 



and $jAf, operators originally introduced by Serot fTHj in his treatment of 1/M^ corrections 
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to currents and charges. (Operators that differ from their standard forms because of sym- 
metrization are denoted here by a tilde. The underlying currents from which these operators 
are derived have a simple behavior under exchange of bra and ket that, unfortunately, is some- 
times lost in the conventional definitions of the weak nuclear operators. The re-definitions 
used here simply restore that symmetry, which is important in dividing elastic operators 
into those that either break or preserve time-reversal invariance.) The single-nucleon matrix 
elements of these operators are real, with the conventions of [18], and transform as follows 

{n{l l/2)j\\dj\\n'{l' 1/2)/) = {-ly'-^ {n'{l' l/2)/||0/|n(/ l/2)j), 6 = M, A, S', S", I-', $" 

(98) 
where n, I, s = 1/2, and j are the nodal quantum number, orbital angular momentum, 
spin, and total angular momentum of the nucleon state. This interchange property is related 
to the time-reversal character of the operator, as we discuss below. A Mathematica script, 
generalized from [25], to evaluate harmonic oscillator matrix elements of these and all other 
operators discussed in this appendix is available from the authors. 



D.l Matrix Elements and the One-body Density Matrix 

Here we address the evaluation of the nuclear matrix elements, including their isospin de- 
pendence and the restrictions imposed by time reversal and parity. We represent operators 
as ajOj, where aj is the overall coupling strength. All nuclear matrix elements appearing in 
the final results of section 3 are singly reduced - reduced in angular momentum only. Thus 
making the suppressed isospin now explicit, and assuming that the nuclear ground state has 
good isospin (thus ignoring charge symmetry breaking) 

A A 

(J,;T,MtJ| J]a,Oj(^) || Jf^T^Mr,) = {J^T^Mt, \\ ^(af=° + af=V3(2)) 6^(2) || J^^TiMr,) 

i=l i=l 

„T=0 ^ T=l Af ^ 

where :: denotes a many-body matrix element reduced in both angular momentum and isospin. 
If the operator Oj is a one-body operator, then 

A 

{J,-T, ■■:■. Y^d'jii) :: J^-Ti) = J^ vl/^f =°(|a|, |/3|) (|a| \\dj\\ 



V2j2^i:i'='M,mM\\o. 



Hm 



{J^T, ;; Y.O'j{i)r{i) ■■ JuT,) = Y. vl/^f=^(|a|,|/3|) (|a| llO.r;; 
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J;T=1 



(|a|,|/3|) (|a| ||6j| 



(100) 



where {\a\ \\0j\\ \/3\) and {\a\ •.-.Ojt:: |/3|) represent single-particle space/spin matrix elements 
reduced in angular momentum and space/spin/isospin matrix elements reduced in both an- 
gular momentum and isospin. The one-body density matrix coefficients "$ r'.- (|a|, |/3|) are 
labeled by the transition i ^ f, angular momentum J and isospin T carried by the transition 
operator, and by the nonmagnetic single-particle quantum numbers \a\ and |/3|. (That is, 
(3 = {n,l,j,mj} for a space/spin matrix element or {n,l,j,mj,mt} for a space/spin/isospin 
matrix element, while |/3| = {n, l,j}.) The sums extend over a complete set of single-particle 



matrix elements. Equation (100) is exact for a one-body operator, regardless of the com- 



plexity of the many-body wave functions. It factors the many-body matrix element of such 
an operator into a sum single-particle matrix elements multiplying numerical coefficients, the 
density matrix. All of the operator physics is contained in the former, while all of the nuclear 
structure physics is isolated in the latter. 



The singly- reduced matrix elements (a;||0j||/3) needed for the present dark matter elastic 
scattering can be conveniently evaluated for a harmonic oscillator single-particle basis, where 
1/3) = \n{l l/2)jmj). The results can be evaluated analytically, producing form factors of the 
form e~^p{y) where y = (g6/2)^, b is the oscillator parameter, and p{y) is a polynomial. 



D.1.1 Symmetry constraints 

The density matrix is defined in second-quantization as 



,J:T 



*;;ai«i 



[jm 



Jf-Tf-: 



«H ® m 



J J\T 



■■ 'Jij ^ i) 



(101) 



where a\. 



\,mj,mt 



('_l)i^-'nj+l/2-mi( 



'mj,~mt 



. With the conventions used here [TTl [18] the 



density matrix is real for parity- and time-reversal-symmetric interactions. Consequently, for 
elastic transitions it follows that 



,J;T 



,J-^T, 



^•r(i«u/3i) = ^-r(i/3i,i«i)(-i)^'^-^'' 



(102) 



so that Eq. (100) can be rewritten as 

A 



[Ji;T,\\ Y,d'j{t)::J,-T,) 



i=l 



v^E^S^^°(i"i'i/5|)o((i«iiio.i 



H,l/3| 



+ (-l)J"--?/3 



IOjII |a|) 



[J,-T,-:: J]6^(0r(^);; J,;7^) 
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v^ E *;r='(i«i' i/5i) ^ ((i«i 110.11 m + i-iy-^'m \\oj\\ i«i)) (103) 



HM 



Now all of the operators as conventionally defined for semi-leptonic weak interactions pTl [TS| 
TU] have real matrix elements with [T71 dS] phase conventions. The contributing operators 
used here have been defined so that they operate symmetrically on bra and ket. Consequently 
it can be shown that the single-particle matrix elements transform as [18] 



(|«| IIOjII |/3|) = (-1)^(1/31 IIOjII |«|) (104) 

where A is either ja + jfs or ja — jp. This gives us the immediate result that the only 
electroweak operators that can contribute to elastic transitions are those transforming with 
A = Jq, — ip, which determines the subset of operators introduced at the beginning of this 
appendix: Mj, Aj, Ej, E'j, $j, and $". If we consider the constraints of parity separately, 
this set is further restricted to the even multipoles of Mj, $j, and $" and the odd multipoles 
of Sj, Sj, and <l>j. This end result does not depend on conventions. 



D.1.2 Nuclear structure 



While in principle Eq. (100) is an exact expression for the many-body matrix elements 
of the contributing operators, in practice, for the complex nuclei of interest to dark matter 
experiments, the density matrix must be taken from a nuclear model. This entails a truncation 
of the sums over |a| and |/3| to some finite Hilbert space, e.g., the active valence shells of a shell 
model, as well as the use of an effective interaction within that valence space, to correct the 
bare nuclear Hamiltonian for the effects of the truncation. The field lacks exact techniques 
for calculating the effective interaction or the corresponding corrections to the operators: 
in practice, most work involves phenomenological adjustments. Some qualitative checks on 
the resulting density matrix can be made by calculating observables other than dark matter 
responses, such as electron scattering elastic charge and (if Jj > 1/2) magnetic form factors 
(or, in the long wavelength limit, the nuclear magnetic moment). 

In the present study, the density matrices needed for the targets we considered came from 
shell-model calculations employing moderately large bases and effective interactions that were 
either determined from global fits to data or based on theoretical input, such as realistic g- 
matrices, supplemented by empirically terms, adjusted to experiment. Potential tests of 
these density matrices include observables such as charge and magnetic elastic form factors 
and magnetic and quadrupole moments. 



D.1.3 Non-recoil Contributions 

The effective theory immediately identifies the Galilean-invariant operators, which in the 
point nucleus limit define all responses connected with target recoil. As explained in section 
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3, this leaves contributions associated with vj^, that is, connected with the intrinsic degrees of 
freedom within the nucleus that would normally be described by introducing the (A-1) Jacobi 
coordinates. Under the assumption that the dark matter couples to to nucleus through local 
one-body interactions with the nucleons within the nucleus, the desired non-recoil intrinsic 
operators can generally be derived. A specific example was given in section 3, where the 
recoil and intrinsic components of the axial charge operator were explicitly identify. Similar 
identifications are possible for all of the interactions where this separation should be done. 

In fact, common practice is to avoid such a separation because the explicit intrinsic op- 
erators tend to be far more complicated than their original one-body forms. Instead, the 
separation is addressed through the nuclear wave functions. While the wave functions are 
also evaluated in an over-complete basis consisting of single-particle Slater determinants, nu- 
merical means are available to remove the center-of-mass motion in the harmonic oscillator 
shell model. This is done by projecting out the center-of-mass motion by adding a large 
multiple of the center-of-mass Hamiltonian to the interaction. The resulting wave functions 
then have the center-of-mass in a Is state: additional "copies" of intrinsic states where the 
center-of-mass is in an excited state are "blown out." This technique requires that calcula- 
tion be done in separable bases: some of the calculations reported here were done in such 
bases, but others required basis truncations that destroy the separability. Furthermore, even 
if the separation is done exactly, the use of the one-body operator forms between such states 
does not provide the exact non-recoil result. Although center-of-mass components in the 
one-body operator cannot generate excitations, they do generate an unwanted center-of-mass 
form factor. 

In fact, all of these considerations are rather technical given other uncertainties in nuclear 
wave functions. The identification of the Galilean-invariant recoil corrections - a trivial task 
in the effective theory - is in most applications more important than whether the nuclear 
physics calculations are done with operators that are precisely intrinsic. The discussion of 
section 3 provided an illustrative example: the entire axial-charge contribution to the elastic 
scattering of dark matter was due to vj^: there are no contributions from vj^ due to the 
combined constraints of time reversal invariance and parity. 

D.2 Operators Needed for Inelastic DM Scattering 

Certain target nuclei may have thresholds for exciting low-lying states that are sufficiently 
low that dark matter interactions can be inelastic. The full scattering probability, analogous 



to Eq. (41 ), will be presented elsewhere. The initial and final nuclear states are then distinct, 
symmetry under time reversal no longer constrains operators or their multi-polarities (though 
parity selection rules remain, depending on the definite parity of the final state, and are used in 
some simplifications below). Here we point out, for those familiar with standard formulations 
of semi-leptonic weak interactions, that a larger set of multipole operators must be defined. 



52 



including the remaining three of the seven basic electroweak operators of [T71 [TH] 

Ajm(?^1 = -z{^VxMf,(gf)|.lv 

SjM(gx) = Mfj{qx)-a 

njMiqx) = njMiqx) + ^SV(g^) = MjMiqx)a " "V + ^ -^MjMiqx) \ ■ a. (105) 

(Note that the symmetrized operator Q is identical to Q' defined and tabulated in [TH]. We 
introduce the alternative notation to avoid confusion, as Q' is defined in [12] as another 
operator. Here we consistently designate with a tilde all operators that require symmetrization 
to restore a simple behavior of matrix elements under interchange of bra and ket.) In addition, 
two symmetrized operators related to Serot's ^jm and A'jj^j are also needed 

^jM{qx) = <^jM{qx) - I^mAqS) = iMfj{qx) ' i^ ^\^) + l]^^ ^ Mfj{qx)\ ■ a 

K"jj,MS) = AjaAq^) - ^MjMiqx) = I -VMjMiqx) | ■ -V - ^Mjuiqx) (106) 

The single-nucleon matrix elements of these operators are again real, under the conventions 
of [IB], but transform with a different sign from those needed for elastic DM scattering. 



{n{l l/2)j\\dj\\n'{l' 1/2)/) = {-iy'^^n'{l' l/2)f\\dj\\n{l l/2)j), 6 = A', S, n, $, A". 

(107) 
This sign is related to operator behavior under time reversal, and thus is the reason that these 
operators do not appear in our treatment of elastic DM interactions. 
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